A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

07(08.01.3; 08.05.1; 08.06.3; 02.12.1; 02.12.3; 02.20.1) 



ASTRONOMY 

AND 
ASTROPHYSICS 



o 
o 
o 



(N ■ 
> 

oo . 
oo , 

^ ■ 

(N ■ 

o : 
o . 

o ■ 

i!^ ■ 

Oh 

6 



Fundamental parameters of Galactic luminous OB stars 
V. The effect of microturbulence* 

M.R. Villamariz^ and A. Herrero^'^ 

^ Institute de Astrofi'sica de Canarias, E-38200 La Laguna, Tenerife, Spain 

^ Departamento de Astrofi'sica, Universidad de La Laguna, Avda. Astroffsico Francisco Sanchez, s/n, E-38071 La Laguna, Spain 



15 December 1999 / 14 March 2000 



Abstract. We study the effect of microturbulence in the 
line formation calculations of H and He lines, in the pa- 
rameter range typical for O and early B stars. We are spe- 
cially interested in its effect on the determination of stellar 
parameters: Tes, logg and specially on the He abundance. 

We first analyze the behaviour of H and He model 
lines between 4 000 and 5 000 A with microturbulence 
and find that for O stars only He i lines and He ii A 4686 
are sensibly affected by microturbulence, and that models 
with lower gravities, the ones suitable for supergiants, are 
more sensitive to it. 

Using a test procedure we show that the expected 
changes in Toff, \ogg and Helium abundance due to the 
inclusion of microturbulence in the analysis, are small. 
We analyze five stars (two late, one intermediate and two 
early O stars) using microturbulence velocities of and 
15 kms~^ and confirm the result of the previous test. The 
parameters obtained for 15 kms~^ differ from the ones at 
kms~^ within the limits of the standard error box of 
our analysis. Only later types reduce their He abundance, 
by 0.02 in e. Comparing with values in the literature we 
find that the range of our changes agree with previous 
results. In some cases other effects can add to microtur- 
bulence, and further reduce the He abundance up to 0.04. 
The quality of the line fits only improves for Hei A 4471, 
but not to the extent of completely solving the so-called 
dilution effect. 

Therefore our conclusion is that microturbulence is af- 
fecting the derivation of stellar parameters, but its effect is 
comparable to the adopted uncertainties. Thus it can re- 
duce moderate He overabundances and solve line fit qual- 
ity differences, but it cannot explain by itself large He 
overabundances in O stars. 



Key words: Stars: atmospheres - Stars: early-type - 
Stars: fundamental parameters - Line: formation - Line: 
profiles - Turbulence 



1. Introduction 

Although microturbulence is a fundamental parameter 
when deriving abundances from stellar spectra, it has been 
systematically ignored when deriving He abundances from 
quantitative spectroscopical NLTE analyses of OB stars 
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(Herrero et al. |1992| , Smith & Howarth |1994| ), bei ng the 
first exceptions to this the works by McErlean et al. (1998) 
and Smith & Howarth ( |199^ ). 

One of the main reasons to ignore it is that the use 
of NLTE techniques reduced the need of microturbulence 
for the reproduction of the observed metallic line strengths 
and they even made abundances less sensitive to the value 
adopted for microturbulence than in LTE analyses (Becker 
& Butler 119891) . 

In quantitative NLTE analyses of OB stars, the lines 
of H and He are used to determine stellar parameters 
(namely, effective temperature, Tcs, logarithmic surface 
gravity, log 5, and He abundance). Their line profiles are 
dominated by the Stark broadening, and microturbulence, 
included in the standard way, only adds an extra Doppler 
width to the thermal broadening. First values for micro- 
turbulence found in LTE for main seq uence B stars, of 
the order of 5 kms~^ (Hardorp & Scholz 1970 ), showed to 
be small enough to be of no importance when compared 
to H and He thermal velocities, that in such hot atmo- 
spheres are well above 5 kms~^. This together with the 
use of NLTE made the influence of microturbulence in the 
overall profile seem negligible. 

When measuring microturbulence in OB supergiants 
the situation is quite different, as its value has always 
showed to be comparable or well above the thermal veloci- 
ties of H and He, and even, in some cases, above the speed 
of soun d in these atmospheres (Lamers 1972| , Lennon & 
Dufton 1986 ). This was not changed b y late r NLTE anal- 
yses (Lennon et al. 1991 , Hubeny et al. 1991, Gies & Lam- 
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bert 1992, Smartt et al. 1997), in which the values of mi 



croturbulence obtained are usuaUy reduced but never to 
values clearly without contradiction. Kudritzki ( 1992| ) and 
Lamers & Achmad (1994) explained that this could be due 
to the presence of wind outflow in these stars, that can 
mimic the effect of microturbulence in the line profiles. 

Most of the previously referred works are on early B 
and late O supergiants. Little work has been done on early 
O stars (see Hubeny et al. 1991), mainly because for them, 
with scarce metallic lines, it is very difhcult to measure the 
value of microturbulence. And it is in the whole range of 
O stars that we are interested. 

One problem that is systematically found in all analy- 
ses of OB spectra is the imposibility of finding a consistent 
fit of all He I lines with a unique set of parameters. This 
difficulty appears worse in supergiants than in main se- 
quence stars, and it is bigger between results from singlet 
and from triplet lines, but also appears within different 
lines in one system. We call this the Hei lines problem. 

This discrepancy between singlet and triplet lines of 
He I was investigated by Voels et al. ( 1989 ), who consid- 
ered that atmospheric extension was responsible and so 
they called this effect "generalized dilution effect". But 
recent works with spherical, mass losing models show that 
this can not be the only reason (Herrero et al., 1995, 
2000|). Furthermore, McErlean et al. (|1998D and Smith & 



abundances during the H-burning phase, this mixing be- 
ing induced by different physical mechanisms like rota- 



Howarth (1998) find that the inclusion of microturbulence 
in the line formation calculations reduces, although not to 
all its extent, the discrepancy between the fits of singlet 
and triplet lines. They also find better and more consistent 
fits for the whole set of Hei lines when microturbulence 
is considered. 

One more reason that encourages us to study the effect 
of microturbulence is that a careful inspection of several 
works shows that even for main sequence stars, values of 
microturbulence comparable to the thermal velocities of 
at least He are obtained in NLTE (Gies & Lambert 1992, 



Kilian et al. 1991), which can invalidate the hypothesis of 



negligible influence of microturbulence in the profiles. 

So in this paper we want to study microturbulence 
in the range of parameters typical for O stars of any lu- 
minosity class, specifically pointing to its effect on the 
determination of stellar parameters and also on the He i 
lines problem. We are specially interested in investigating 
whether the inclusion of microturbulence can solve the 



He discrepancy (Herrero et al. 1992), as Vrancken (1997), 
McErlean et al. ( |1998| ) and Smith & Howarth ( |1998D sug- 
gest. This now well known problem has induced a lot of 
theoretical improvements in both evolutionary and model 
atmosphere theories in order to explain it. Recent works 
including sphericity and mass loss in the spectral analysis 



(Herrero et al. 1995, 200C, Israehan et al. 200C) show that 



with these new models, more suitable for these stars than 
the plane-parallel hydrostatic ones, the discrepancy is 
not solved. Evolutionary models including additional mix- 
ing inside the star can explain enhanced He photospheric 



Maeder, 1997, Maeder & Zahn, 



tion and turbulent diffusion (Denniscnkov, 1994, Meynet 

|1998|) or turbulent dif- 



fusion and semiconvection (Langer 1992| ). Sometimes an- 



gular rotational velocity changes during evolution are in- 

Heger, 



eluded (Langer & Heger, 1998 



1998). For a recent 



review, see Maeder & Meynet (200C). Therefore it becomes 
of great importance to see whether microturbulence can 
be responsible for the He discrepancy. 

In Sect. 2 we present line formation calculations of H 
and He lines and their behaviour with microturbulence in 
the O stars domain. In Sect. 3 we determine the effect of 
including microturbulence in the determination of stellar 
parameters, and in Sect. 4 we analyze some stars with 
and without considering microturbulence. Sections 5 and 
6 are then dedicated to our discussion and conclusions, 
respectively . 

2. Microturbulence in H and He line formation 
calculations 

In order to study the behaviour of H and He lines with 
microturbulence, we perform line formation calculations 
in the parameter range typical for O to early B stars of 
any luminosity class: 

-between 30 000 and 45 000 K in Teg 
-between 3.05 and 4.00 in logg (in c.g.s. units) 
-for e= 0.10 and 0.20, with e = jv(j^)^^(g) 
where N(x) is the number density of atom X. We fol- 
low the classical technique of calculating a NLTE model 
atmosphere of H and He, in radiative and hydrostatic equi- 
librium and wit h plan e-parallel geometry (calculated with 
ALI, see Kunze, 1995), and then solve the statistical equi- 
librium and transfer equations and perform the formal so- 
lution for the lines of H and H e cho sen using DETAIL & 
SURFACE (Butler & Giddings |l985| ). In this final step we 
have included UV metal line opacities in the calculations, 
in order to obtain more r ealistic profiles (see Herrero, 1994 , 
and Herrero et al., 2000 , for details). As shown in the last 
reference, plane-parallel models are unable to reproduce 
properly spectra of massive OB stars around 50 000 K and 
hotter, that is why we stop our study at 45 000K. 

Microturbulence is introduced in the standard way, by 
adding an extra Doppler width to the thermal broadening 
of the line, which is then convolved with the rest of the 
broadening mechanisms. We have considered it in both the 
equations of statistical equilibrium and radiative transfer. 
As we don't consider other motions (turbulent or not) in 
the determination of the stellar structure, we prefer to re- 
strict the introduction of microturbulence to the absorp- 
tion coefficient. Therefore, we do not include it in the cal- 
culation of the structure of the atmosphere via turbulent 
pressure. For similar reasons also we do not even bother 
about its dependence on depth. We don't treat separately 
the effect on the populations and line profiles, as this has 
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been recently studied by McErlean et al. (1998) and Smith 
& Howarth ( |1998D (whose results we support). 

We make line formation calculations for the lines we 
usually consider to perform our analyses: Hj and for 
Hi; He I AA 4387, 4922, 4471 A for Hei; and Heii AA 
4200, 4541, 4686 A for Hen. 

In order to study the behaviour of these lines with 
microturbulence, we perform line formation calculations 
for microturbulence velocity values from to 20 kms~^, 
for different sets of parameters, representing the spectral 
types and luminosity classes we are interested in: 

- Teff = 30 000, log 5= 3.05, e= 0.10 for late O and early 
B supergiants. The same with \ogg= 4.00 for dwarfs. 

- Teff= 35 000, \ogg= 3.20, e= 0.10 for "middle" O 
supergiants. The same with log 5= 4.00 for dwarfs. 

- reft= 40 000, logg= 3.40, e= 0.10 for early O super- 
giants. The same with logg= 4.00 for dwarfs. 

- Teff= 45 000, log 5= 3.60, 0.10 for very early O 
supergiants. The same with log 5= 4.00 for dwarfs. 

The values of log g for supergiants are close to the low- 
est that could be converged for every T^s. For dwarfs of 
increasing Tcff, log 5 becomes a little below 4.00, but this 
value still represents well this luminosity class. Finally, 
these calculations were repeated for e= 0.20, to look for 
differential effects with He abundance. 

Figs. 0, H and ^ display the behaviour of H and He line 
profiles with microturbulence, for three different models, 
adequate respectively for late O supergiants, late O dwarfs 
and early O supergiants. In Table |l| equivalent widths are 
given for all these lines. The first two Figs, will allow us to 
compare the effects of microturbulence when varying log 5, 
whereas the first and last Figs, will be used to compare 
the effects when varying TeS . The other calculated models 
behave consistently with what it is shown in Figs. ^ ^ and 
will not be further discussed. 

Looking first at Fig. |l| we see that there are two differ- 
ent line groups. The first one is composed by the strong 
H lines and the relatively weak lines Hell AA 4200, 4541. 
All these lines are not affected by microturbulence, al- 
though our highest values are comparable to the thermal 
velocity of H atoms and well above that of He atoms. 
The reason is that the Stark broadening dominates the 
profiles, and therefore masks the changes that microtur- 
bulence produce. 

The second group is composed by Hei lines and the 
strong He 11 A4686 line. He i lines show sensitivity to micro- 
turbulence in both core and wings, while the core of He 11 
A4686 is desaturated by microturbulence, but its effect is 
masked in the wings by the Stark broadening, because, 
as it has been explained by Smith & Howarth (1998), the 



effect of microturbulence will depend on the stepness of 
the line wings. 

When we then compare Figs. and |2| we see that the 
increased pressure broadening (directly related to the in- 
creased density) simply reduces the effect of microturbu- 
lence in the line profiles. Thus, as gravity increases, the 




4470 4471 4472 4473 
Wovelength [Angstrom) 



4685.0 4635.5 4686.0 4686.5 4587.0 
Wavelength (Angstnonn) 



Fig. 1. From top to bottom and left to right: , Hp , Hei 
AA 4387, 4922, Hen AA 4200, 4541, Hei A 4471, Heii A 
4686, for the model with Tes= 30 000, log g= 3.05, e= 0.10, 
representing late O and early B supergiants. Microturbu- 
lence takes values of (solid fine), 5 (dotted), 10 (dashed), 
15 (dash dotted) and 20 (dash tri-dotted) kms""'^. 



effect of microturbulence decreases for all considered lines. 
In the case of Hell lines this effect is reinforced by the 
displacement introduced in the He ionization equilibrium, 
that makes He ii lines weaker. 

More interesting is the comparison with higher tem- 
peratures. Looking to Fig. |l| and Fig. || we can see the 
behaviour of the lines of supergiants with increasing T^g. 
For higher values of T^s marginal effects on H i and He II 
line cores are seen (Fig. but only Hei lines and Hell 
A4686 are again really sensitive to microturbulence. He ii 
line profiles are still insensitive, although they are now 
much stronger and even have a larger equivalent width 
than Hei lines, because they are dominated by the Stark 
broadening. The reason why He I line profiles are still sen- 
sitive to microturbulence, although they are much weaker, 
has to be found in the infiuence of microturbulence on the 
shape of the absorption profile of weak lines. Thus, we see 
that Hei AA 4387, 4922 are not saturated in the whole 
temperature range. Their equivalent widths increase with 
increasing microturbulence, until, as Toff increases, they 
become weak enough to have equivalent widths, but not 
line profiles, independent of microturbulence. This hap- 
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Table 1. Equivalent widths of the lines in Figs. |l], || and ^, in mA. For every line the values for microturbulence 
velocities of 0, 5, 10, 15 and 20 kms~^ are given from top to bottom. Model parameters in the first column are Tcff in 
K, logg in cgs units and e. 



Params. 




H,3 


He I 4387 


He I 4922 


He I 4471 


He II 4200 


Hen 4541 


Hen 4686 


30 000 


1895 


2149 


403 


503 


622 


195 


234 


416 


3.05 


1896 


2150 


415 


526 


637 


195 


234 


418 


0.10 


1896 


2151 


446 


585 


682 


196 


235 


415 




1904 


2160 


485 


664 


750 


199 


238 


410 




1904 


2173 


525 


753 


836 


201 


241 


401 


30 000 


3831 


3832 


599 


733 


1158 


42 


50 


161 


4.00 


3832 


3833 


605 


747 


1167 


42 


50 


160 


0.10 


3834 


3836 


618 


780 


1189 


42 


50 


157 




3839 


3843 


636 


825 


1220 


42 


50 


151 




3850 


3857 


654 


873 


1252 


42 


50 


145 


40 000 


2005 


2393 


120 


238 


391 


491 


614 


685 


3.40 


2004 


2392 


118 


240 


402 


491 


613 


684 


0.10 


2001 


2385 


118 


255 


434 


492 


613 


685 




2020 


2406 


120 


279 


480 


519 


646 


667 




2011 


2392 


120 


299 


530 


522 


649 


669 




pens at 40 000 K for Hei A 4387 (see Table |l|) and at 
45 000 K for Hei A 4922. For them, Doppler broadening is 
important enough to let microturbulence shape the pro- 
files appreciably in the whole temperature range. Hei A 
4471 is saturated at T^s 30000 K, but it desaturates at 
Teff= 35 000, and starts behaving like the other Hei lines. 



Finally, we point out that all these results also apply 
for models with e = 0.20, and that they are a natural 
extension to O stars of the results of Mc Erlean et al. 
(|1998D and Smith & Howarth (|199^). 
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3. The influence of microturbulence in the 
determination of parameters of O-stars 

In this work we are interested in estimating the changes 
introduced in the derived stellar parameters by the inclu- 
sion of microturbulence in the analysis of O stars, with 
special attention to its influence in the derived He abun- 
dance. As we saw in the previous section, models of low 
gravity, corresponding to supergiants, are more sensitive 
to microturbulence, and so we expect supergiants to show 
the largest changes. 

Values of microturbulence of 10, 12 kms~^ are found 
for late O and early B stars (see references in the intro- 
duction), therefore we decided to choose a fixed value for 
microturbulence of 15 kms~^ to perform a test of the in- 
fluence of microturbulence in the determination of Toff , 
log (7 and He abundance. 

As a first approach to the problem, we take a model 
spectrum at ^=15 kms~^ as the observed spectrum, to 
which we fit model profiles in a grid around it at ^=0 
kms~^, with parameters differing in 500 to 1000 K in Tcs, 
0.05 to 0.10 in logg and 0.02 to 0.04 in e (these smah 
values are suggested by previous inspection of a larger, 
coarser grid). To maximize the effect of microturbulence, 
we consider neither rotational broadening nor instrumen- 
tal broadening. Taking advantage of the fact that we are 
using only theoretical profiles (although the one with ^= 
15 kms""'^ has been adopted as the observation) we use a 
least-squares fit procedure to determine the best fit. 

For each line of each model in the grid, we calculate 
the quadratic difference with the observed line, and then 
we add the results for all lines of a given ion. For example 
for H we calculate: 



Table 2. xj^ of models of the grid that best fit the observed 



xf(H) 



y - 

H lines 



X,obs } 



as the result for H for model i of the grid, where wx is the 
spectral sampling and AA is the wavelength interval con- 
taining the whole line. We perform these calculations for 
all the models in the grid and then normalize these values 
to their minimum, so that a 1.00 gives us the model that 
best fits the lines of the ion at all temperatures, gravities 
and He abundances. The calculation is performed in the 
same way for He lines, taking Hei and Heii lines sepa- 
rately. Finally we define for model i: 
Xf= x?(H)+xf(Hei)+xf(Heii) 

so that the best fit is adopted to be that of the model 
with the minimum value of X\. 

We have carried on the exercise for a model with 
r^ff=40 000 K, logg=3.40, e=0.10 and ^=15 kms'^ as the 
observed spectrum. In Fig. || we have plotted contour lev- 
els of in the T^a - logg plane for e = 0.10. We see 
that, at this He abundance, there are two models that 
fit well the observation, with parameters Tcff^ 40 000K, 
\ogg= 3.40 and T^s^ 42 000K, \ogg= 3.45. Although the 
last model has a slightly lower value of xf (5.190), the 



spectrum 


with parameters 


reff=40 000 


K, logg: 


=3.40, 


e=0.10 and ^= 


15 kms 


^, see 


text for details. 




Tcfs 


logs 


e 


xim 


xmei) 


xmell) 


xi 


41000 


3.50 


0.06 


9.46 


2.50 


9.12 


21.07 


41000 


3.45 


0.08 


1.51 


2.25 


1.09 


4.85 


40 000 


3.40 


0.10 


1.36 


1.40 


2.67 


5.44 


42 000 


3.45 


0.10 


1.00 


3.19 


1.00 


5.19 


39 000 


3.35 


0.12 


2.07 


1.82 


2.06 


5.96 


40 000 


3.35 


0.14 


6.81 


1.60 


7.13 


15.54 




TEFF/IOOO 



Fig. 4. X Contour levels in the Tcff - logg plane for e= 
0.10, showing that the best fit for the observed model at 
this abundance is still a model with the same parameters 



dispersion in the individual xf values is larger, refiecting 
the fact that looking at the line fits we would choose the 
other as the best compromise. This is a indication of the 
differences between both criteria, we should improve our 
X^ fitting criterium in order to match the fitting by eye. 

In Tabic ^ wc list the results for the best fits at other 
He abundances. Wc sec inmcdiatcly that for e — 0.06 and 
0.14 we find large values of x^, indicating a relatively poor 
fit. On the contrary, the models for e= 0.12, 0.10, 0.08 fit 
with approximately the same quality (without a more de- 
tailed study we cannot say whether the difference in x^ is 
significant). In a fit by eye we would choose the model at 
e 0.08 as the best fit. Actually, this is the model with the 
lowest value of a-nd has a lower He abundance than the 
observed spectrum. The reason is that the larger gravity 
rather than the He abundance mimics the microturbulence 
effect. In principle, it would be also possible to choose the 
model with e= 0.12, in which the higher He abundance 
corresponds with a lower gravity. Fig. || shows the fit of 
the "best model" to the observed spectrum. The lack of ro- 
tation and instrumental broadening makes the differences 
in He I lines appear evident. 

Thus, the conclusion of our exercise is that neglection 
of microturbulence will change slightly the derived stellar 
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0,90 
> O.BO 
' 0.70 
- 0.60 




To analyse stellar spectra we first determine the pro- 
jected rotational velocity of the star (see Herrero et al. 



4334 4335 4338 4340 



4B56 4838 4860 4862 4854 4866 




4470 4471 4472 4473 
Wavelenglh (Angstrom) 



■474 4684.5 4685.0 4685.5 4586.0 4686.5 4687.0 
Wovelsngth (Angstrom) 



Fig. 5. The fit of the "best model" (dotted, Tcs= 41 000 
log 5= 3.45 e= 0.08, ^= kms^^ ) to the observed one 
(solid, Teff= 40 000 logg= 3.40 e= 0.10, ^= 15 kms'^ ). 



parameters, but keeping them within our standard error 
box (see next Sect.). The direction in which parameters 
are moved will depend on the criteria used for defining a 
model as "the best model fit" , such as the relative weight 
given to different lines, or to the core and wings of a line. 
Thus the definition of these criteria for the "best model" 
will be a critical point of any future automatic fitting pro- 
cedure, which will soon be demanded by new observing 
capabilities. 



4. Spectral analysis including microturbulence 

Now we want to see the difference we obtain in the de- 
rived parameters of real observed O stars spectra when we 
include microturbulence in the calculations. We selected 
two late-0 supergiants to see how our results match with 
previous works, as well as one intermediate and two early 
ones, in order to cover the whole O spectral type. We chose 
HD 5 689 in particular, an 06 V fast rotator with a high 
He-overabundance derived in prev ious w orks without mi- 
croturbulence (see Herrero et al., 2000), to check how it 
can modify this strong overabundance. 

The description of the obser vations and data reduction 
can be found in Herrero et al. (|l992| , |200C|) . 



1992 for details). This is an additional parameter that, in 
fact, reduces the effect of microturbulence in the deriva- 
tion of stellar parameters. The results of the analyses are 
listed in Table |. 

The fitting procedure and adopted errors of ± 1 500 
K in Toff, ± 0.1 in logg and ± 0.03 in e are explained in 
Herrero et al. ( |1999[ ). 

To make the new analyses with microturbulence we 
construct a model with the same parameters but for ^= 
15 kms^^, and also a small grid around it with changes 
in Toff of ± 500 to 2000 K, in logg of ± 0.05 to 0.15 
dex and in e of ± 0.02 to 0.04. These small changes are 
suggested by the results of the preceeding section, and are 
in fact confirmed by the small differences in the fit we 
see between the two models at and 15 kms~^. The new 
best fit is found as explained above, and it gives the new 
parameters for the star at ^=15 kms""'^. 

All the stars except HDE 338 926 have been previously 
analyzed by our group without microturbulence (Herrero 
et al., 1992, 200C). In the present analysis, small differ- 
ences for the stars in common, specially with respect to 
the first reference, can be found due to the larger weight 
given here to Hei A4387, rather than to Hei A4922, and 
to the inclusion of line-blocking in the calculations (the 
temperatures of the late type supergiants are slightly hot- 
ter in the present analysis). All this also helps to solve 
the small difference between He i A 4387 and He i A 4922 



found in Herrero et al. (1992) for these late O supergiants. 

In Table I we see the results for the five stars. Changes 
in the parameters induced by microturbulence are not be- 
yond the standard error box of the analysis , and in par- 
ticular, we see that He abundance is reduced for four of 
the stars, but only slightly. In Table ^ we see the values 
obtained for the radius, mass and luminosity, following the 
same procedures outlined in Herrero et al. ( |1992D . We see 
that changes are not significant. (We point out that the 
values for mass, radius and luminosity given in Table I for 
HD 210 809 and HP 18 4 09 differ slightly from those given 
in Herrero et al. (1992). These differences are not due to 
the new parameters, but to a change in the distanc e mod - 
uli. Herrero et al. used the values from Humphreys ( 1978 ), 
while we use the values from Garmany & Stencel ( 1992|) ). 

In Figs. ^ to ^ we can see the fits with microturbu- 
lence and 15 kms^^^. The quality of the fits to individual 
lines is the same at both values, although the triplet line 
He I A 4471 becomes stronger at 15 kms^^ , which produces 
an improved fit, but not to the extent as to completely re- 
solve the so-called dilution effect. Sometimes, a compro- 
mise between Hei A4387 and Hei A4922 is adopted, like 
in HD5 689 and HD 210 809 (realize that for 15 kms"! 
the first one is slightly weak and the second one slightly 
strong) . 

With respect to HDE 338 926, that has been analyzed 
here for the first time, we have to point out that the anal- 
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Table 3. Stellar parameters with values of microturbulence of and 15 kms^^. Effective temperatures are given in K 
and surface gravities in c.g.s. units. Spectral classifications are taken from Walborn (1973), except that of HD5 689, 
which is taken from Garmany & Vacca (1991), and that of HDE338 926, which is taken from Humphreys ( |1978| 
Humpreys lists this star as BD+24 3866). Some final parameters are extrapolated as we could not converge the model 
with the indicated values of logg but others with \ogg 0.05 larger, with all other parameters the same. This is indicated 
by giving the parameters in italics 



Star 


Clasif. 


Vr sin i 






log 


',9 




e 






kms 


kms"^ 


15 kms^^ 


kms"^ 


15 kms^^ 


kms"^ 


15 kms^^ 


HD 14 947 


05 If+ 


140 


45 000 


45000 


3.50 


3.45 


0.15 


0.15 


HD 5 689 


06 V 


250 


40 000 


40 000 


3.40 


3.35 


0.25 


0.22 


HDE 338 926 


08 f 


130 


34000 


34500 


3.00 


3.00 


0.15 


0.13 


HD 210 809 


09 lab 


120 


33 500 


34 500 


3.10 


3.15 


0.10 


0.08 


HD 18 409 


09.7 lb 


160 


31500 


32 000 


3.10 


3.10 


0.11 


0.09 



Table 4. Radius, mass and luminosity for the analyzed objects. As in Table y, italic numbers mean values obtained 
by extrapolation 



Star 


kms 


R/Rq 

15 kms"^ 


kms" 


M/Mq 

15 kms~^ 


log(L/L0) 
kms"^ 15 kms-^ 


HD 14 947 


14.8 


14.9 


26.5 


24.1 


5.91 


5.91 


HD 5 689 


7.7 


7.8 


7.8 


7.4 


5.13 


5.14 


HDE 338 926 


24.4 


24.3 


23.8 


23.5 


5.86 


5.88 


HD 210 809 


18.2 


17.9 


16.0 


17.4 


5.57 


5.61 


HD 18 409 


18.8 


18.8 


18.7 


18.6 


5.50 


5.52 



ysis has been very difficult. The final parameters are ex- 
trapolated beyond the models we could converge. How- 
ever, we are confident that these parameters do charac- 
terize the star appropriately (or as appropriately as those 
of other stars), because the fits with already converged 
models are reasonably good, and because we are able to 
fit the star with converged models when we do not con- 
sider line-blocking. Thus, it is only the small temperature 
increase introduced by line-blocking that moves the star 
beyond the convergence region. The value we obtain for 
the evolutionary mass, derived from the tracks by Schaller 



et al., (1992), is 57.6 Mq without microturbulence, and 
56.8 MfT) with 15 kms~^. Comparing with the values in 
Table W, it is clear that HDE 338 926 also shows a mass 
discrepancy, as usual. Of course, we will have to reana- 
lyze HDE 338 926 in the future with spherical, mass lossing 
models. 

5. Discussion 

We will not go here into the discussion of the real physi- 
cal entity of microturbulence, neither on the validity of the 
aproximation of small scale turbulent movements assumed 
to introduce it just as an extra Doppler width, which may 
be not suitable for the big values we deal with. This is 
beyond the scope of this work. As explained in Sect. 2 
we do not introduce microturbulence in the structure cal- 
culations because we think that, having neglected other 
motions, it is actually not more physically consistent to 
introduce a turbulent pressure term. We just accept the 
necessity of using microturbulence in the analysis of stel- 



lar spectra, specially in the determination of metal abun- 
dances, that is our final interest. 

We find that contrary to our previous considerations. 
He I lines and He 11 A 4686 do have profiles sensitive to the 
usual values of microturbulence found in OB stars. We 
confirm this for the whole range of parameters describing 
O and early B spectra, in agreement with McErlean et 
al. ( |1998| ) and Smith & Howarth ( |1998| ) for early B and 
late O supergiants respectively. This sensitivity is first due 
to the fact that the Stark broadening is not dominating 
completely these profiles and second to the high values 
of microturbulence involved, that are comparable to the 
thermal velocity of He ions. Taking (0.84 T^s) as repre- 
sentative of the temperature in the zone of formation of 
the lines, we find that for roff= 30 000 K wth of He ions 
is 10 kms-\ and for Tes= 45 000 K it is 12 kms'^ For 
H ions, with a fourth of the He atomic weigth, thermal 
velocities are two times larger, so thermal broadening and 
Stark broadening dominate the profile. For all other Hell 
lines Stark broadening is hiding the effect of microturbu- 
lence. 

Quantifying the effect of microturbulence on the deter- 
mination of stellar parameters we find that they are not 
changed beyond the standard error box of our analyses. 
This implies that for stars with high He overabundances, 
such as HD 5 689 analysed here, the He discrepancy will 
not be solved by considering microturbulence. This result 
is also supported by what we obtained in Sect. 3. 

This seems to be in contradiction with Smith & 



Howarth (|1998| ) and McErlean et al. (|1998| ), who afhrm 
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Fig. 6. Fits with ^ = kms ^ (solid line) and ^ = 15 
kms^^ (dotted) to the star HD 14 947. The parameters 
are T^s= 45 000, log.g= 3.50, e= 0.15 at ^ = kms"! and 
the same at 15 km^^, which are the converged models 
closest to the final adopted parameters, with log g = 3.45 



Fig. 7. Fits with ^ = kms ^ (soUd hne) and ^ = 15 
kms^^ (dotted) to the star HD 5 689. The parameters are 
Teff^ 40 000, log.9= 3.40, e= 0.25 at ^ = kms"! and 
reff= 40000, log5= 3.35, e= 0.22 at ^ = 15 kms^^ 



that solar He abundances are found for the supergiants 
they analyse, when microturbulence is considered. How- 
ever, the situation is still unclear, as an analysis of pub- 
lished results can show. 

Let us begin with the He overabundances obtained by 
Herrero et al. (|1992D . Following the argument by Smith 
& Howarth ( |1998[ ) and McErlean et al. ( |1998[ ), the pre- 
ferred use of the strong line He i A 4922 in the analyses 
(slightly more sensitive to microturbulence than He i A 
4387) would be responsible for the high He abundances 
found in that work. As a conclusion, the derived He over- 
abundances would be an artifact introduced by the neglec- 
tion of microturbulence. 

This argument applies in the cases in which Herrero 
et al. ( 1992 ) found a discrepancy between Hei A 4922 and 
He I A 4387, but not in the rest of the cases. 

In Fig. ^ we have plotted He abundance versus Toff, 
as derived by Herrero et al. (1992), using different symbols 



for stars for which a discrepancy between these two lines 
was reported. We can see that the priority given to He i A 



4922 can affect some of the abundances derived, specially 
at lower temperatures. 

To illustrate the situation, let us briefly discuss one of 
these objects for which Herrero et al. report a difference 
between He i A 4922 and He i A 4387: HD 210 809. This 
star has been analyzed here again, where we have given 
more weight to He i A 4387 (contrary to Herrero et al.). 
The final difference between the He abundance obtained 
here with a microturbulence of 15 kms~^ and that ob- 
tained by Herrero et al. is 0.04 (our new He abundance 
being lower). This, however, has to be adscribed to dif- 
ferent effects. First, giving more weight to He I A 4387 in 
absence of microturbulence results in a higher tempera- 
ture (as already stated by Herrero et al.), which leads to 
a lower He abundance to keep the fit of Hen lines. The 
effect of line blocking, not included in Herrero et al., goes 
in the same direction. This accounts for a difference of 
0.02 and the additional difference of 0.02 is purely due to 
microturbulence (see Ta ble ^ ) . Comparing our Fig. ^ with 
Fig. 5 of Herrero et al. ( 1992 ) we see that the effects con- 
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Fig. 8. Fits with ^ = kms ^ (solid line) and ^ = 15 
kms^^ (dotted) to the star HDE 338 926. The parameters 
are reff=34000, log.9= 3.05, e= 0.15 at ^ = kms"! and 
Teff= 34500, log5= 3.10, e= 0.13 at ^ = 15 kms^i. Again 
these are the closest models that could be converged, for 
a final adopted logg of 3.00 



sidered here help to reduce the discrepancy between the 
two He I lines. 

The discrepancy bet ween H e i A 4387 and He i A 4922 
found by Herrero et al. ( |1992[ ) is larger at lower temper- 
atures, and thus, around 30 000 K it cannot be solved 
by varying T^s or introducing line blocking. Therefore, 
around or below this temperature microturbulence is at 
present the only considered effect that could bring them 
into agreement. 

In the literature we can find analyses of early B and 
late O supegiants including microturbulence (Gies & Lam- 
bert^ |1992| , Smith & Howarth, |1994 |l998| , Smith et al., 
1998| , McErlean et al., |199^ ). In Table |^ we hst stars for 



which parameters have been derived with and without mi- 
croturbulence. 

Comparisons in Table |5| have to be done with care. 
It mixes results from different authors and different cri- 
teria. We see however that it indicates that the changes 
found here are consistent with those of other authors, i.e., 



Fig. 9. Fits with ^ = kms ^ (solid line) and ^ = 15 
kms~^ (dotted) to the star HD 210809. The parameters 
are roff=33 500, log.9=3.10, e=0.10 at ^ = kms"! and 
Teff=34 500, log5=3.15, e=0.08 at ^ = 15 kms"! 



stellar parameters are changed within the uncertainties 
adopted here. These changes do not follow a clear, sys- 
tematic pattern (i.e., going always in the same direction 
when introducing microturbulence) and thus we conclude 
that the effect of microturbulence is indeed not larger than 
present-day uncertainties. An exception to this might be 
K Ori, which shows a large change in He abundance. How- 
ever, this large reduction of the He abundance found by 
McErlean et al. ( 1998| ) as compared to Lennon et al. ( 1991 ) 
is accompanied by a large change in the stellar parame- 
ters. Even if we attribute the whole change to the effect 
of microturbulence (and not to a new analysis with new 
data and more refined calculations) we see that k Ori is by 
far the coolest of the stars in Tables || and || and thus we 
expect a larger influence of microturbulence for it. In addi- 



tion, we note that McErlean et al. (199?) do not exclude a 
larger He abundance for this star. We conclude that data 
in the literature do not lead to the conclusion that the He 
discrepancy in O and early B supergiants is completely 
due to microturbulence, although microturbulence helps 
to reduce it. 
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Table 5. OB supergiants whith parameters determined with and without microturbulence. Values of microturbulence 
vary between 10 and 15 kms^^. References are as follows: (1) Herrero et al. ( 2000D ; (2) this work; (3) Herrero et al. 
(|l99|); (4) Smith & Howarth ( |l994| ); (5) Smith et al. (|l99|); (6) Smith & Howarth ( |l99^ ); (7) Lennon et al. ( pMl| ); 
(8) McErlean et al. (|l99|) 



Star 


Clasif. 


Vi sin i 




ff 


log 


9 




e 




Ref. 






kms ~^ 


C = o 




e = o 




C = o 




e = 


e/0 


HD14 947 


05 If+ 


140 


45.0 


45.0 


3.50 


3.45 


0.15 


0.15 


1 


2 


HD 5 689 


06 V 


250 


40.0 


40.0 


3.40 


3.35 


0.25 


0.22 


1 


2 


HD 210 809 


09 lab 


100 


33.0 


34.5 


3.10 


3.15 


0.12 


0.08 


3 


2 


HD 154 368 


09.5 lab 


85 


33.0 


32.0 


3.07 


3.00 


0.13 


0.13 


4 


5 


HD 123 008 


ON9.5 lab 


90 


33.5 


33.0 


3.07 


3.05 


0.17 


0.15 


4 


5 


HD 152 003 


09.7 lab 


80 


30.8 


29.7 


2.90 


3.10 


0.09 


0.09 


6 


6 










30.5 




3.00 




0.09 




5 


HD18 409 


09.7 lb 


160 


31.0 


32.0 


3.10 


3.10 


0.12 


0.09 


3 


2 


HD 154 811 


OC9.7 lab 


125 


31.5 


31.0 


3.15 


3.10 


0.09 


0.09 


4 


5 


K Ori 


BO. 5 la 


80 


25.0 


27.5 


2.70 


3.00 


0.20 


0.10 


7 


8 
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Fig. 10. Fits with ^ = kms ^ (solid line) and ^ = 15 
kms~^ (dotted) to the star HD18 409. The parameters are 
Teff= 31500, logg= 3.10, e= 0.11 at ^ = kms^^ and 
Tee= 32 000, \ogg= 3.10, e= 0.09 at ^ = 15 kms"!. We 
could not use neither H^ nor Hei A 4922 for the fit, as we 
don't have observations of this star in this range. 

Except for the discussed reduction of the He abun- 
dance, we do not see a clear pattern in the parameter 
changes in Tables | and |. Tes and log g can both increase 
or decrease when microturbulence is introduced, and thus 
we are tempted to conclude that these changes are still 
dominated by internal inconsistencies in the analyses that 
appear when we compare values that cannot be distin- 
guished within the adopted uncertainties. 



0,00 L 

45 



EFFECTIVE TEMPERATURE 



Fig. 11. e versus Tes (in thousands of Kelvin) in the sam- 
ple of Herrero et al. ( |1992 ). Filled diamonds represent 
stars for which Herrero et al. report a large discrepancy 
between He i A 4922 and He i A 4387. Open diamonds rep- 
resent stars for which Herrero et al. report a moderate 
discrepancy between these two lines. Plus signs represent 
stars for which Herrero et al. do not report a discrepancy 



About the problem of the consistency of the fits to 
He I lines, we find that the dilution of He i A 4471 is only 
partially solved, even when line-blocking is considered as 
here. The fits to the rest of He I lines do not improve much 
either in the stars we analyze here. So the consideration 
of both microturbulence and line-blocking in the analysis 
cannot completely make an agreement between the results 
from triplet and singlet Hei lines, though they help to 
improve it. 

6. Conclusions 

We study for the first time the effect of microturbulence 
in the whole O spectral range, from late to early O types. 
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Introducing microturbulence in the solution of the sta- 
tistical equilibrium and transfer equations, and then in 
the formal solution (i.e., in the absorption coefficient), we 
conclude that for higher gravities the effect in the lines is 
negligible. This together with the low values of microtur- 
bulence usually found for high luminosity class stars, lead 
us to conclude that only O supergiants have sensitivity to 
microturbulence effects at a given Tes- 

In examining the behaviour with microturbulence of 
the lines we use in our analysis, we show that only Hei 
lines and the core of He ii A 4686 A are sensitive to mi- 
croturbulence. For He I lines we show that, as should be 
expected, there is not a constant pattern for each line, 
but it depends on the parameters considered, that deter- 
mine the strength of the line and its degree of saturation. 
This invalidates generalizations to the whole O spectral 
range made upon results obtained just for a certain spec- 
tral type. 

In order to quantify the sensitivity of stellar parame- 
ters to microturbulence we find that changes in the param- 
eters induced by a value of microturbulence of 15 kms~^ 
are enclosed within the standard error box of our analyses. 
We think that the lack of a clear pattern in the changes 
induced in Toff and log 5 is just due to the fact that we 
are varying the parameters within this error box. We do 
however find a systematic change in e towards lower He 
abundances when microturbulence is introduced. In par- 
ticular we find that late O supergiants show a decrease of 
0.02-0.04 in e (this last value when including other effects 
that add to microturbulence), which is in agreement with 
previous results pointing to the inverse relation between 
the microturbulence assumed and the He abundance ob- 
tained (Smith & Howarth 1998| ). Early types are less sen- 
sitive to microturbulence, and might not show a difference 
in the derived He abundance. 

Thus microturbulence is not capable of explaining the 
He discrepancy at all for early O stars, and neither it is 
for late O types with high overabundances. 

Looking to individual lines we find that the fits to He i 
A 4471 A are improved when considering microturbu- 
lence, but not to the extent of completely explaining its 
dilution. On the other hand Hei AA 4922, 4387 A are 
sometimes slightly better and sometimes slightly worse 
fitted, in the last case with model cores a bit too strong or 
a bit too weak, respectively. The rest of the lines keep the 
same quality in the fit. The Hei lines problem is there- 
fore only partially solved by simply considering microtur- 
bulence, even with line-blocking included in the model 
profiles. 

Therefore our conclusion is that microturbulence is af- 
fecting the derivation of stellar parameters, but its effect is 
comparable to the adopted uncertainties. Thus it can re- 
duce moderate He overabundances and solve line fit qual- 
ity differences, but it cannot explain by itself large He 
overabundances in O stars, and we are forced to conclude 
that these are due to other effects, whether real or caused 



by artifacts in our analyses. This last point will probably 
not find a definitive answer until we are able to derive re- 
liable abundances of C, N and O in the atmospheres of O 
stars that we can correlate with the He abundances. 
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